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fi]  Observations  of  airglow  at  630  nm  (red  line)  and  557.7  nm  (green  line)  during  the 
February  2002  campaign  at  the  High  Frequency  Active  Auroral  Research  Program 
(HAARP)  heating  facility  are  analyzed.  We  find  that  during  injections  toward  magnetic 
zenith  (MZ)  the  green  and  red  lines  gain  ~5  R  within  ~1  s  and  ~20  R  within  ~10  s, 
respectively.  We  term  this  period  the  onset  of  the  HF-induced  airglow.  A  model  of  the 
onset  at  magnetic  zenith  is  developed.  It  accounts  for  background  photoelectrons  and 
dissociative  recombination  of  0\.  It  is  shown  that  heating  and  acceleration  of  background 
electrons  dominate  the  airglow  onset.  We  propose  a  scenario  for  the  generation  of  strong 
Langmuir  turbulence  for  injections  outside  the  Spitze  region,  including  magnetic 
zenith.  INDEX  TERMS:  2403  Ionosphere;  Active  experiments;  2483  Ionosphere:  Wave/particle 
interactions;  2471  Ionosphere:  Plasma  waves  and  instabilities;  2481  Ionosphere:  Topside  ionosphere; 
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1.  Introduction 

[2]  A  distinctive  feature  of  HF  modification  experiments 
is  the  excitation  of  airglow  at  630.0  and  557.7  nm  by  high- 
power,  high-frequency  (HF)  radio  waves  [e.g.,  Sipler  et  al., 
1974;  Bernhardt  et  al.,  1989;  Pedersen  and  Carlson,  2001; 
Gustavsson  et  al.,  2001,  2002;  Kosch  et  al.,  2000,  2002). 
Enhancements  up  to  ~500  R  (Rayleighs)  with  the  green-to- 
red  ratio  C,gr  as  high  as  >0.3  have  been  reported 
[Gustavsson  et  al.,  2002;  Kosch  et  al.,  2002;  Pedersen  et 
al.,  2003).  The  excitation  energies  ex  of  the  O('D)  and 
O('S)  states  responsible  for  the  red  and  green  lines  are  er  = 
1.96  eV  and  eg  =  4.17  eV,  respectively.  The  population  of 
energetic,  e  >  Z\,  electrons  can  increase  significantly  due  to 
stochastic  and  resonant  interactions  with  plasma  turbulence 
generated  by  a  heating  wave  [e.g.,  Gurevich  et  al.,  1985; 
Dimant  et  al.,  1992;  Mantas  and  Carlson,  1996;  Gurevich 
and  Milikh,  1997;  Istomin  and  Leyser,  2003).  The  former 
raises  the  electron  temperature  Te,  the  latter  accelerates  a 
group  of  electrons  in  the  high-energy  tail  of  the  initial 
distribution.  Both  effects  are  well  documented  [e.g., 
Carlson  et  al,  1982;  Gustavsson  et  al.,  2001). 

[3]  Gustavsson  et  al.  [2002]  emphasized  that  to  interpret 
C^r  >  0.1  in  terms  of  electron  heating  requires  unrealizable 
Te>  2  eV,  pointing  out  the  importance  of  electron  acceler¬ 
ation.  It  is  commonly  believed  that  electrons  are  most 
efficiently  accelerated  by  Langmuir  ( l)  turbulence.  The 
generation  of  Langmuir  waves  is  usually  described  in  terms 
of  nonlinear  instabilities  of  ordinary  (o)  mode  pump  waves, 
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known  as  the  parametric  decay  (PDI)  or  oscillating  two- 
stream  (OTSI)  instabilities  [Fejer,  1979).  In  the  presence  of 
background  suprathermal  (e  »  Te)  electrons,  many  more 
energetic  electrons  are  accelerated  than  would  be  in  Max¬ 
wellian  plasmas  [ Mishin  and  Telegin,  1986). 

[4]  Recent  observations  show  that  the  airglow  maximizes 
during  HF  injections  toward  magnetic  zenith  [Kosch  et  al., 
2002;  Pedersen  et  al.,  2003).  The  same  is  true  for  the 
intensity  of  Langmuir  waves  [Isham  et  al.,  1999]  and 
electron  heating  [Rietveld  et  al.,  2003]  observed  by  the 
EISCAT  UHF  radar.  Furthermore,  the  red  line  is  excited  at 
magnetic  zenith  even  at  extremely  low  effective  radiative 
power  P0  ~  2  MW  [Pedersen  et  al.,  2003]. 

[5]  Ordinary  mode  waves  with  incident  angles  0  outside 
the  Spitze  region,  0  >  0C,  reflect  at  altitudes  Ha  below  the 
standard  reflection  altitude  H0  where  the  local  plasma 
frequency  fp  ~  10 4v/«7  Hz  equals  the  driver  frequency /0. 

Here  0C  =  arcsin(^^sinx),  fc  is  the  local  electron  gyro- 

frequency,  ne  is  the  electron  density  in  cm-3,  and  x  is  the 
magnetic  dip  angle  [e.g.,  Mjolhus,  1990].  At  HAARP  x  — 
14.5°  and ~  1.4  MHz  at  altitudes  near  200  km,  so  that  for 
fo  =  l  MHz  the  Spitze  angle  is  0<-  ~  5.9°. 

[6]  Figure  1  shows  a  schematic  of  ray  trajectories  for 
ordinary  HF  waves  injected  vertically  and  toward  magnetic 
zenith  into  a  horizontally  stratified  ionosphere  at  HAARP. 
Obliquely  incident  radiation  does  not  form  standing-wave 
patterns  and  swelling  is  absent.  Thus,  given  P0  =  150  MW 
and  distance  R  =  250  km,  the  wave  amplitude  is  E0  ~ 
4.7 -JP^/R  ~  0.2  Vim.  With  /0  =  7  MHz  and  Te  =  0.1 2eV, 
the  HF  energy  density  at  the  reflection  point  is  W0  =  ^  ~ 
10 ~*neTe,  sufficient  to  drive  the  PDI/OTSI  at  H0  [Fejer, 
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Figure  1.  Schematic  of  ray  propagation  for  ordinary 
HF  waves  injected  toward  local  vertical  and  magnetic 
zenith.  The  magnetic  field  B  direction  is  indicated  by 
arrows.  The  heights  of  reflection  and  upper  hybrid 
resonance  are  shown  by  horizontal  lines.  A  light  dashed 
line  shows  the  Spitze  angle  direction.  The  regions  of 
excitation  of  Langmuir  turbulence  appear  as  bold  lines  near 
the  reflection  points. 


1979].  However,  mismatch  of  frequencies  at  Ha  <  H0  suppress 
these  instabilities  for  injections  outside  the  Spitze  region. 

[7]  Gurevich  et  al.  [2002]  suggested  that  decreased 
plasma  densities  within  striations  generated  by  heating 
waves  permit  necessary  phase  matching  [cf.  Muldrew, 
1978].  Changes  in  the  refraction  index  (self-focusing)  due 
to  striations  develop  within  tens  of  seconds  after  turn-on  and 
explain  some  features  of  the  spatial  distribution  of  the  HF- 
induced  airglow.  However,  the  rise  time  of  Langmuir  waves 
is  ~10  ms  [Isham  et  al.,  1999].  Besides,  striations  are 
generated  in  the  upper  hybrid  layer  [e.g.,  Vaskov  et  al, 
1981;  Lee  and  Kuo,  1983].  Hence  the  reflection  height  of  the 
heating  wave  must  be  at  or  above  this  height.  For  magnetic- 
zenith  injections  at  HAARP,  this  can  be  satisfied  only  if 
/o  <  5.4  MHz.  However,  the  strongest  airglow  at  HAARP 
occurred  with  higher  jo  values  [Pedersen  et  al.,  2003]. 

[8]  Kuo  et  al.  [1997]  showed  that  Langmuir  waves  can  be 
excited  by  upper  hybrid  ( uh )  waves  with  amplitudes  exceed¬ 
ing  the  threshold  value  E'uh  ~  0.15  V/m.  Such  uh- waves  can 
be  generated  through  the  linear  conversion  of  the  o-mode  on 
pre-existing  field-aligned  irregularities  [  Wongetal. ,  1981]  or 
through  the  parametric  decay  o  — >  uh  +  Ih  [e.g.,  Istomin  and 
Leyser,  1995];  Ih  stands  for  lower  hybrid  waves.  Linear 
conversion  within  the  upper  hybrid  layer  proceeds  with  no 
set  threshold.  However,  the  parametric  decay  develops  if  E0  > 
Eq1'  ~  I.6/0  mV/m,  provided  that  the  inequality  0.015/ 
\fQ  —  s  ■  fc |  <  0.5  MHz  (integer  s  >  3)  is  valid.  Here  f0 
stands  for  the  heating  frequency  in  MHz.  The  rise  time  of 
the  uh- wave  is  Tuh  ~  1-3  ms.  Consistent  with  observations 
[Isham  et  al.,  1999],  this  makes  the  generation  of  Langmuir 
turbulence  possible  within  ~10  ms. 

[9]  Besides  electron  impact,  dissociative  recombination 
(DR)  of  0\  is  an  efficient  source  of  the  0('D)  and  0(lS) 
states.  Since  the  rate  of  dissociative  recombination  decreases 
with  Te,  it  is  usually  not  considered  in  the  theory  of  the  HF- 
induced  airglow.  However,  the  quantum  yield  for  the  O('S) 


state  from  dissociative  recombination  grows  significantly 
with  Te  and  with  the  vibrational  temperature  of  oxygen  ions 
[Guberman,  1997;  Peverall  et  al.,  2001],  Furthermore, 
charge  exchange  0+  +  02  — >  02  +  O  is  the  dominant  source 
of  O2  production  in  the  nighttime  F  region.  Its  rate  increases 
significantly  whenever  02  is  vibrationally  excited  [Viggiano 
and  Williams,  2001],  We  expect  a  high  degree  of  excitation 
of  molecular  species  in  the  HF-illuminated  region  and  thus 
enhanced  yields  of  the  green  line  emissions. 

[10]  A  consistent  theory  of  the  HF-induced  airglow 
accounting  for  the  ion/neutral  chemistry,  modification  of 
the  heated  spot,  and  self-focusing  has  yet  to  be  worked  out. 
We  emphasize  that  responsible  processes  develop  within  tens 
of  seconds.  However,  considerable  airglow  appears  within 
first  few  seconds  after  the  HF  transmitter  turns  on.  We 
designate  this  interval  as  the  onset  of  HF-induced  airglow. 

[11]  This  paper  develops  a  model  for  the  onset  of 
airglow  enhancements  at  magnetic  zenith  accounting  for 
the  roles  of  0\  dissociative  recombination  and  energetic 
photoelectrons.  The  following  section  describes  the  onset 
characteristics  from  the  HF  heating  experiments  at 
HAARP.  Section  3  describes  the  model.  In  particular,  we 
propose  a  scenario  (section  3.3)  for  the  generation  of 
strong  Langmuir  turbulence  outside  the  Spitze  region 
including  magnetic  zenith  with  upper  hybrid  waves  as 
the  primary  source.  The  final  section  compares  modeling 
results  with  optical  measurements. 


2.  Airglow  Onset:  HAARP,  February  2002 

[12]  The  HAARP  facility  is  located  outside  of  Gakona 
Alaska  (62.4°N,  145.15°W).  During  the  course  of  HF 
heating  experiments  between  03:00  and  05:00  UT  in  Feb¬ 
ruary  2002  several  passes  of  Defense  Meteorological  Sat¬ 
ellite  Program  (DMSP)  satellites  flew  to  the  east  and  west  of 
the  HAARP  location.  Each  of  these  spacecraft  carries  a  pair 
of  upward  looking  particle  spectrometers  designed  to  mea¬ 
sure  fluxes  of  downcoming  electrons  and  ions  with  energies 
between  30  eV  and  30  keV.  Consistent  with  prevailing  quiet 
geomagnetic  (Kp  =  2)  conditions,  the  equatorward  boundary 
of  auroral  precipitation  was  several  degrees  in  latitude 
poleward  of  Gakona.  However,  both  before  and  after  the 
local  time  of  sunset,  the  spectrometers  detected  fluxes 
of  Rdowncoming  electrons  whose  spectra  monotonically 
decreased  with  energies  between  30  and  100  eV.  These 
subauroral  fluxes  consist  of  photoelectrons  that  originated 
in  the  still  sunlit  southern  ionosphere. 

[13]  Intense  green-line  emissions  were  observed  between 
03:48  and  05:00  UT  (17:48-19:00  LT)  on  February  13, 
2002  events  of  the  HAARP  optics  campaign.  Readers  are 
referred  to  the  report  of  Pedersen  et  al.  [2003]  for  graphic 
examples  of  HAARP-  induced  airglow  and  the  intensities  of 
red/green-line  emissions  recorded  during  the  period  of 
interest.  O-mode  waves  were  injected  toward  magnetic 
zenith  at  f0  =  7.8  MHz  and  at  full  power  of  0.94  MW 
( Pq  —  1 65  MW).  The  transmitter  was  programmed  to  turn  on 
for  5  minutes  exactly  on  the  minute.  Subsequent  pulses 
followed  5  minutes  pauses.  However,  the  fourth  pulse  in  the 
sequence  started  at  04:21:18  UT.  A  transmitter  problem 
produced  a  false  start,  04:30:00-04:31:00  UT,  prior  to  the 
long  pulse  that  began  at  04:32:12  UT.  In  the  course  of  this 
long  pulse,  the  critical  frequency  of  the  F  layer  dropped 
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Figure  2.  Heater-induced  green  (top  panel)  and  red 
(bottom  panel)  lines  during  February  13,  2002  injections. 
Vertical  solid  and  dash-dotted  lines  indicate  heater  turn-ons 
and  turn-offs,  respectively. 
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below  the  transmitter  frequency  and  the  sequence  was 
terminated. 

[14]  To  account  for  natural  (background)  variations,  a 
polynomial  fit  was  made  to  intensities  measured  along  the 
magnetic  meridian  in  each  all-sky  image,  with  the  heated 
area  blocked  out.  Throughout  the  studied  period  the  differ¬ 
ence  between  the  background  airglow  b\  and  that  of  the 
heated  volume  4>  during  heater-off  periods,  was  less  than 
~3-5  and  ~7 — 10  R  for  green-  and  red-line  emissions, 
respectively.  This  small  difference  determines  the  accuracy 
of  heater-induced  airglow  measurements  A \  =  I\  —  b\. 

[15]  Figure  2  presents  a  detailed  view  of  the  last  three  HF- 
induced  enhancements  of  February  13,  2002  at  f0  = 
7.8  MHz.  Data  points  are  from  two  red  line  exposures 
compiled  at  0  and  24  s  after  the  beginning  of  each  minute, 
and  green-line  exposure  beginning  at  12  s  into  the  minute. 
After  04:00  UT,  the  exposure  duration  was  7.5  s.  By 
examining  cases  when  the  heater  turned  on  part  way 
through  an  exposure,  we  see  that  the  green-line  emissions 
increased  by  ~5  (~10)  R  within  ~1.5  (~7.5)  s,  while  the 
red  line  gained  25  (  30)  R  within  ~13.5  (~20.5)  s. 

[16]  Starting  at  04:27:15  UT  on  February  9,  2002,  a  P0  ~ 
124  MW  wave  at  a  frequency  6.8  MHz  was  injected  toward 
magnetic  zenith.  Figure  3  shows  red-  and  green-line  emis¬ 
sions  before  and  during  the  pulse  turn  on  from  the  heated 


volume  and  the  background  region.  These  data  show  that 
the  difference  between  natural  green-line  emissions  from 
magnetic  zenith  and  ~70  km  apart  was  less  than  1  R.  The 
same  is  true  for  the  red-line  emission  before  the  heater 
turned  on.  This  suggests  that  the  background  airglow  was 
nearly  uniform.  When  the  heater  turned  on,  green  (red)-line 
intensities  increased  by  ~5  (~30)  R  within  ~4.5  (~16.5)  s. 
It  is  relevant  to  note  that  Kosch  et  al.  [2002]  reported 
~500  R  enhancements  of  the  red-  and  green-line  emissions 
produced  by  the  EISCAT  superheater  (P0  ~  600  MW) 
within  the  first  5-s  frame. 

3.  Modeling  Airglow  Onset 

[17]  For  modeling  purposes  we  assume  that  the  major 
species  constituents  of  the  neutral  atmosphere  and  iono¬ 
spheric  plasma  are  well  described  by  the  MSIS-E  [Hedin, 
1991]  and  IRI  [ Bilitza ,  2001]  models,  respectively.  Although 
photodissociation  of  02  accounts  for  ~50%  of  the  back¬ 
ground  airglow,  it  is  unaffected  by  heating  and  hence  is  not 
included.  Basic  processes  and  their  rate  coefficients  are 
discussed  at  length  by  Solomon  et  al.  [1988],  Rees  [1989], 
Witasse  et  al.  [1999],  and  Fox  and  Sung  [2001],  Unless 
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Figure  3.  Heater-induced  green  (top)  and  red  (bottom) 
lines  during  February  9,  2002  injections.  Dots  and  asterisks 
stand  for  emissions  from  the  magnetic  zenith  and  back¬ 
ground  regions,  respectively.  Vertical  dashed  line  indicates 
the  heater  turn-on  time. 
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otherwise  noted,  we  use  the  rate  coefficients  of  Fox  and 
Sung  [2001], 

[is]  As  mentioned  above,  radiative  emissions  from  a 
heated  volume  reflect  the  superposed  results  of  complex 
wave,  plasma,  and  chemical  interactions.  To  render  our 
modeling  of  this  complex  chain  of  interactions  intelligible, 
we  divide  this  section  into  three  parts.  The  first  subsection 
concerns  the  basic  contributors  to  the  natural  and  artificially 
excited  radiation  budgets.  The  subsection  on  the  back¬ 
ground  ionosphere  estimates  the  distributions  of  ion  species 
and  energetic  photoelectrons  present  at  heater  altitudes  prior 
to  turn  on.  The  third  subsection  describes  a  flow  chart 
encapsulating  our  concept  of  how  heater-injected  energy  is 
transformed  into  the  various  wave  modes  that  heat  and/or 
accelerate  ambient  electrons  that  interact  with  ambient 
neutrals  to  produce  onset  green-  and  red-line  emissions. 

3.1.  Basic  Processes 

[19]  The  red  (r)  and  green  (g)  line-photons  are  emitted  by 
atomic  oxygen  in  the  transition  from  the  0(}D)  and  0(lS) 
states  to  0(3P )  and  0(lD)  states,  respectively.  In  photo¬ 
chemical  equilibrium,  the  volume  emission  rate  r|x  is 
calculated  from 


Figure  4.  Cross-sections  of  the  0(lD)  and  0('S)  excita¬ 
tion  and  the  Schumann-Runge  dissociation  by  direct 
electron  impact.  Squares  and  dots  are  the  measured  cross- 
sections  from  Doering  and  Gulcicek  [1989];  solid  lines  are 
interpolations  by  Majeed  and  Strickland  [1997]. 


<■ 


■  A\  ■  [Ox]  =A\ 


Q\ 


L\  + 


(1) 


Here  [Ox]  stands  for  the  density  of  0('D)  or  OfS)  in  cm-3; 
Q\  and  L\  are  the  corresponding  production  and  loss  rates, 
respectively;  Ar  ~  0.007,  A&  ~  0.009,  Ag  ~  1 .2,  and  A&  ~ 
1.3  are  the  Einstein  transition  probabilities  in  s_I.  The 
column  emission  rate  (in  R)  is  found  by  integrating  (l) 
along  the  line  of  sight 


4iHx  =  10~6  [*'  r\fz)dz  (2) 

J’n 


[20]  The  rates  of  electron  impact  excitation  O  +  e  — > 
O0D,'S)  +  e  and  the  Schumann-Runge  dissociation  02  + 
e  — >  O  +  0(lD)  +  e  are  calculated  as  follows 

2x  =  [°1  ’  *x  =  \°\  ■  4lt  /  <7X(e)$(e)fife  (3) 

Here  $(e)  is  the  electron  differential  number  flux,  xx,  ex, 
and  ax  are  the  rate  coefficients,  threshold  energies,  and 
cross-sections,  respectively.  We  employ  electron  impact 
cross-sections  suggested  by  Majeed  and  Strickland  [1997]. 
Figure  4  shows  three  of  the  major  electron  impact  cross 
sections.  It  is  worth  noting  that  xer  (Tej  and  x|  (Te)  calculated 
with  a  Maxwellian  distribution  of  thermal  electrons  are 
close  to  the  approximations  used  by  Mantas  and  Carlson 
[1996]  and  Gurevich  and  Milikh  [1997],  respectively. 

[21]  The  dissociative  recombination  of  02  produces 
Q\R  =  K\Rnc[02\.  The  rate  coefficients  are  usually  approx¬ 
imated  as  %\R  =  0x*  ’  O',  where  0x*  are  the  quantum  yields. 
Estimated  yields  for  0(!6)  were  subject  to  substantial 
discrepancy  in  the  literature  for  many  years.  Recent 
experiments  [Kella  et  al.,  1997;  Peverall  et  al. ,  2001; 
A.  Petrignani  et  al.,  Vibrationally  resolved  rate  coefficients 
and  branching  in  the  dissociative  recombination  of  02, 
submitted  to  Journal  of  Chemical  Physics,  2003]  and  theory 


[Guberman,  1988,  1997]  appear  to  have  reconciled  this 
controversy.  They  established  that  a  and  Px*  depend  not 
only  on  Te  but  also  on  the  vibrational  population  n°l  - 
[02(j)]l[02]  where  j  =  0,  1,. . .  designates  2the 


vibrational  level. 

[22]  For  the  ground  vibrational  state,  0^O)  ~  1.15,  a(0)  ~  2 
10-7t7°'63  [cf.  Fox  and  Sung,  2001],  and 


-v)) 


(4) 


for  1  <  Te  <  10,  where  re  =  re[K]/300.  Importantly,  the 
inequalities  |3^'2)  >  (3^0)  and  $1,2)  <  (3[.0)  hold.  For  the 
Maxwell-Boltzmann  distribution  of  N^l  with  vibrational 
temperatures  Tv  ~  0.25  eV,  one  obtains^*’  ~  1.2|3i0)  and 
O.850j.o). 

[23]  Collisional  quenching  0#  +  X  -*  O  +  X  is  important 
for  the  0(lD)  state,  while  for  0('S)  at  altitudes  >150  km  it 
can  be  ignored.  Here  X  stands  for  O,  02,  N2,  or  e,h.  The 
corresponding  loss  rates  are  L*  =  10“ 11  xqx\X],  where  xS  — 
0.65  t°'14,  xS2~  3.2  exp  (0.225/t„),  y.Nq  ~  1 .8  exp  (0.36/t„), 
and  x?  ~  28.7  t®’91.  Figure  5  shows  the  altitude  profile  of 
the  total  loss  rate  calculated  with  the  MSIS-E  and  IRI 
parameters  at  04:30  UT,  February  13,  2002. 

[24]  Finally,  the  radiative  cascade  (RC)  0(lS)  — ►  0(lD)  + 
hv5 577  yields  QRC  =  r^. 

3.2.  Background  Ionosphere 

[25]  Simultaneous  observations  from  a  digisonde  located 
at  the  HAARP  site  were  used  to  correct  the  IRI  model  and 
to  determine  that  the  reflection  height  at  magnetic  zenith 
decreased  from  ~250  to  ~235  km  between  04:00  and 
04:30  UT  on  February  13,  2002.  At  the  same  time,  the 
shadow  height  increased  from  ~1 50  to  ~230  km  indicating 
the  presence  of  photoelectrons  [e.g.,  Doering  et  al.,  1975]. 
Furthermore,  the  southern-hemisphere  region  magnetically 
conjugate  to  HAARP  remained  in  sunlight,  even  after  local 
sunset.  Thus  a  large  fraction  of  the  photoelectron  spectrum 
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0(1  D)  loss  rate,  s  1 
Figure  5.  0(lD)  loss  rate. 


generated  in  the  conjugate  ionosphere  [e.g.,  Peterson  et  al., 
1977]  had  access  to  the  ionosphere  above  HA  ARP.  Electron 
spectrometers  on  DMSP  satellites  observed  the  high-energy 
tail  of  the  conjugate  photoelectron  spectrum. 

[26]  Figure  6  shows  altitude  profiles  of  the  background 
electron  temperature  T0  and  electron  and  0+  densities 
calculated  from  the  corrected  IRI  model.  The  twilight  0\ 
density  above  ~200  km  is  defined  by  the  charge  exchange 
and  dissociative  recombination 


[0+]  ~xc 


[02]  [Q+l 


Here  the  rate  of  charge  exchange  %ce  after  Viggiano  and 
Williams  [2001]  is 


Figure  7.  First,  upper  hybrid  waves  are  generated  due  to 
the  parametric  decay  o  — ■>  uh  +  Ih.  Its  threshold  field  is 
Eoh  ~  0.08  V/m  if  the  matching  conditions  are  met  at 
xuh  =  kuhre  <  1,  where  re  is  the  thermal  electron  gyrora- 
dius  and  kuh  is  the  w/i-wave  vector.  Landau  damping  of 
short-scale  lower  hybrid  waves  raises  the  threshold  other¬ 
wise  [ Mishin  et  al.,  1997].  For  f,  =  7.8  (6.8)  MHz  and  0  ~ 
14.5°,  the  matching  conditions  at  the  reflection  height  H0 
are  easily  satisfied  for  xuh  ~  0.3  (~0. 1 5).  Finally,  the  free 
space  field  of  the  incident  wave  Eo  ~  0.25  (0.2)  V/m 
exceeds  Eoh. 

[29]  When  the  amplitude  of  the  primary  w/i-wave  Euh 

exceeds  ~10  mV/m,  it  excites  other,  lower  frequency, 
w/i-wave  with  the  growth  rate  of  order  yuh  ~  r~hl  [Zhou 
et  al.,  1994],  The  same  is  true  for  subsequently  generated 
w/i-waves.  The  frequency-step  of  this  spectral  transfer  is 
quite  small  ~  ii.\xuhbkuhlkuh  ~  0 .2yjme/mi,  where 

mjmi  represents  the  electron/ion  mass  ratio.  As  many 
spectral  steps  occur  before  the  parametric  instability  satu¬ 
rates,  the  resulting  w/i-energy  spectrum  consists  of  a  large 
number  AS;  10  of  spectral  peaks,  each  of  order  E^/(yuhrruh) 
(well-known  weak  turbulence  “cascading”  [e.g.,  Sagdeev 
et  al.,  1991]).  Thus  the  total  w/i-energy  density  can  be 
estimated  as  W*/A  ~  W0  ~  7.5  (5)  eV/cm3  with  the  r.m.s 
amplitude  E*  ~  ^4-nW*  VaEq. 

[30]  Upper  hybrid  waves  with  amplitudes  Euh  exceeding 
the  threshold  value  E'uh  ^  0  .1  !x]!h  V/m  generate  short- 
scale,  k,  >  3~l/2/re,  Langmuir  waves  [Kuo  et  al,  1997]. 
Due  to  a  weak  dependence  of  the  growth  rate  7,  on  kuh, 
one  can  take  Euh  ~  E *.  As  this  value  greatly  exceeds  the 
threshold  eL  the  growth  rate  can  be  evaluated  as  y/  ~ 
veEuh!Eluh  [Kuo  et  al.,  1997],  which  amounts  to  yt  ~ 
vey/\02xuf,AEQ-  Here  ve  is  the  electron  elastic  (transport) 
collision  frequency. 

[31]  The  growth  of  Langmuir  waves  at  TJTt  <  4  saturates 
via  induced  scattering  by  ions  that  transfers  energy  toward 


r-ce  ^  xj?  (hq’  +  n^[4.9 lg(T/n  +  1)  -  2.5]) 

(-^)).  <6) 

a  is  the  recombination  rate  in  cm3/s, 

[7]  stands  for  the  density  of  species  Y  in  cm-  ,  and  t  = 
1.5t,  +  t„  —  1  with  t iin  =  7]  „[K]/300. 

[27]  The  IRI  model  predicts  Te=  T0  <  0.2  eV  (Figure  6). 
At  these  temperatures  and  ne>  10s  cm-3,  a  Maxwellian 
distribution  FM  is  a  good  representation  of  the  thermal 
electron  distribution  function  (Fth)  in  the  F  region  [Mishin 
et  al.,  2000].  The  differential  number  flux  $5(e)  of  photo¬ 
electrons  in  the  F  region  [e.g.,  Rees,  1989]  at  e  >  zs  -  ~mv2 
~  5  eV  can  be  approximated  by  a  power-law  function 


*,(e)  = 


ps  —  0.5 


where  ns  <  100  cm 


to3  104  ioi3  106 

number  density,  cm 


3.3.  HF-Perturbed  Ionosphere  Figure  6.  Altitude  profiles  of  the  electron  temperature 

[28]  Our  scenario  for  the  excitation  of  plasma  turbulence  (solid  lines),  electron  density  (dashed  line),  and  0+  density 

and  subsequent  electron  energization  in  the  HF-illuminated  (dotted  line).  T0,  Tly  and  T2  stand  for  0,  1,  and  ~1,5  s  after 

region  at  magnetic  zenith  is  represented  schematically  in  turn-on. 
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Figure  7.  Block-diagram  showing  synergy  of  a  model  of 
HF-induced  airglow  at  magnetic  zenith.  Shaded  boxes 
indicate  the  steps  to  be  worked  out. 


small  k  [e.g.,  Zakharov  et  al.,  1976].  This  process  is 
governed  by  the  pump  u/i-wave  power.  In  particular,  if  the 
width  A k)  of  the  parametrically  unstable  /-wave  spectrum 
exceeds  A kc  ~  k,vJh,  the  Langmuir  wave  energy  W, 
accrues  in  a  state  with  k  — ♦  0  (Langmuir  condensate) 
[Zakharov  et  al,  1976;  Zakharov,  1984],  The  spectral  width 
can  be  evaluated  as  Ah  ~  Hfo/fc)2W*/(neTere).  Thus,  for 
E0  >  0.15(/ove)1/6/\/A  £  0.2  V/m  the  “condensate” 
condition  A kt  >  k/vJh  is  fulfilled.  The  dynamics  of  the 
Langmuir  condensate  is  defined  by  the  modulational  insta¬ 
bility  and  collapse  leading  to  establishment  of -strong 
(cavitating)  Langmuir  turbulence  [e.g.,  Galeev  et  al., 
1977;  Zakharov,  1972,  1984], 

[32]  We  emphasize  that  this  applies  not  only  at  the 
reflection  layer  Ho  but  also  well  below  it,  wherever  the 
matching  conditions  for  the  parametric  decay  o  — >  uh  +  Ih 
are  met  and  local  values  of  xuh  <  1.  Given  the  plasma 
density  profile  shown  in  Figure  6  and  f0  ~  6-8  MHz,  the 
altitude  extent  of  this  region  is  ~10  km.  One  should 
compare  this  value  with  a  few  100-m  size  of  the  Airy 


pattern  at  H0,  where  the  PDI/OTSI  instabilities  develop 
during  injections  within  the  Spitze  region  [. Fejer ,  1979],  We 
believe  that  this  is  the  key  factor  in  exciting  strong  airglow 
at  magnetic  zenith. 

3.3.1.  Heating 

[33]  Coliisional  damping  of  high-frequency  o-,  uh-,  and 
/-waves  is  the  major  source  of  (stochastic)  electron  heating. 
Mishin  et  al.  [2000]  showed  that  at  [JV2]  >  107/o  cm-3  the 
process  of  A^-vibrational  excitation  dominates  the  forma¬ 
tion  of  the  distribution  function  of  ionospheric  electrons 
with  energies  3.5  >  e  >  ev»  —  1-8  eV,  which  can  be 
represented  as  follows 


Flh(e  >  Evil, )  - 


'£  +  Eytt 


5K.(£yft) 

<s) 


x  exp 


L  JTTt 


\  /e°-5k(s 

eSifci 

f  S“1/2k(0 di 

Jtvib 


(8) 


Here  k,2(e)  =  v,/(e)/uee(e)  +  0 .5e/7),  where  vjvee  is  the  ratio 
between  electron  inelastic  and  electron-electron  collision 
frequencies.  A  total  HF  wave  energy  density  WHF  -C 
10 ~2n0Te  is  assumed.  At  energies  e  <  £„,■*  the  distribution  is 
close  to  while  at  e  >  3.5  eV  it  may  slightly  deviate  from 
Fm  due  to  the  0(lD)  and  O('S)  excitation. 

[34]  The  electron  temperature  can  be  evaluated  from  the 
energy  balance 


1.5 nedTjdt  =  Te-  148  •  neTe  -  (9) 

where  qe\\  is  the  (parallel)  electron  heat  flux,  Ft.  ~  veWHF  is 
the  volume  heating  rate,  and  6  =  (v,//ve)  is  the  coefficient  of 
inelastic  losses  averaged  over  the  total  distribution  F,h. 

[35]  Given  ve  ~  500  s_1,  is  £40  (25)  keV/cm3/s.  Note 
that  to  match  their  optical  observations,  Mantas  and 
Carlson  [1996]  used  the  values  of  Te  between  ~50  and 
~125  keV/cm3/s  at  ~260  km  for  ^  ~  5  MHz.  Gustavsson 
et  al.  [2001]  used  Tc  ~  60  keV/cm3/s  for^o  ~  4  MHz  to  fit 
the  electron  temperature  profile  maximum  4000  K  at 
~220  km  observed  by  the  EISCAT  UHF  radar.  Neither  of 
them  accounted  for  the  decrease  of  5(7))  due  to  the  deviation 
of  the  thermal  electron  distribution  from  Maxwellian.  Taking 
that  into  account,  the  calculation  results  of  Mantas  and 
Carlson  [1996]  and  Gustavsson  et  al.  [2001]  can  be  scaled 
to  the  Ti  ^-profiles  in  Figure  6,  pertaining  to  our  case. 

3.3.2.  Acceleration 

[36]  Resonant  Ih-  and  /-wave-particle  interaction  accel¬ 
erates  electrons.  Musher  et  al.  [1978,  1986]  analyzed  in 
detail  the  dynamics  of  Ih- waves  excited  by  an  external 
source.  At  1  <  TJTt  <  4,  high-frequency  waves,  fh  > 
fc\J  Intel  mu  are  dominated  by  induced  scattering  by  ions. 
The  rate  of  spectral  transfer  toward  smaller  frequencies  is 
1th  ~  ( fcZ/fhW,hln0Te ,  where  W,h  is  the  energy  density  of 
lower  hybrid  waves.  Equating  the  growth  rate  yuh  to  ylh 
yields  the  energy  density  at  the  saturated  state  WfhlnQTe  ~ 
luhfihlfc2  -  10-5-5. 

[37]  The  dynamics  of  Ih- waves  in  the  low-frequency,  fh  < 
fcy/lmjmu  region  is  dominated  by  the  lower  hybrid 
collapse  [e.g.,  Musher  et  al.,  1978].  The  threshold  energy 
density  is  quite  low  WfhlnTe  ~  xlh(me/m,)(fc/f0)2  <  10~6, 
and  is  surely  exceeded  in  our  case.  In  the  course  of  collapse, 
the  longitudinal  and  transverse  dimensions  of  a  Ih- cavity, 
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/|l  ~  l±y/mi/xu,me,  diminish.  Finally,  the  Ih- wave  energy  in 
a  collapsing  cavity  is  absorbed  by  particles  via  Landau 
damping  [Musher  et  al.,  1978,  1986;  Sotnikov  et  al.,  1978; 
Shapiro  et  al.,  1993].  The  density  of  accelerated  electrons  in 
the  magnetic  field-aligned  high-energy  tail  of  the  distribu¬ 
tion  function,  determined  by  assuming  resonant  particles 
carry  away  all  of  the  energy  pumped  into  collapsing  cavities 
[e.g.,  Shapiro  et  al.,  1993],  is  of  order  <10 ~sne.  This 
//j-contribution  appears  insufficient  to  account  for  the 
observed  HF-induced  airglow.  It  is  worth  noting  that 
//i-collapse  also  leads  to  transversely  accelerated  ions  that 
absorb  the  collapsing  energy  at  approximately  the  same  rate 
as  electrons  [Sotnikov  et  al.,  1978]. 

[38]  Strong  Langmuir  turbulence  consists  of  an  ensemble 
of  collapsing  cavitons  that  transfer  the  wave  energy  toward 
small  scales.  Small-scale,  k  >  iVvmin,  waves  are  absorbed 
by  fast  v  >  vmjn  3>  \J Te /me  electrons  via  Landau  damping, 
thereby  increasing  the  population  of  suprathermal  electrons. 
The  distribution  of  accelerated  electrons  Fa(e)  =  jr$a(z)  in 
weakly  magnetized  fc<.fP  plasmas  can  be  found  from  the 
kinetic  equation  [e.g.,  Galeev  et  al.,  1977] 


&±  =  }_<L\A.  fX  dk 

dt  v2  dv  mnv  J^/v  dv 


(10) 


where  w*/(2ir)  ~’fp  and  Wk  are,  respectively,  the  frequency 
and  spectral  energy  density  of  Langmuir  waves.  The  latter  is 
determined  from  the  requirement  for  energy  balance 

m  d 

dt  +  dk 

Here 

r*--(v)  $>•(¥)  <”> 

is  the  rate  of  Landau  damping;  k(t)  ~  (tc  —  t)~2  is  defined 
by  the  collapse  law  in  the  absorption  region  [ Pelleiter , 
1982],  and  tc  is  the  collapse  time. 

[39]  The  steady  state  solution  of  (10)— (12)  at  e  >  emin  = 
5»ievminis  a  power-law  function  [Pelleiter,  1982;  Shapiro 
and  Shevchenko,  1984]  that  yields  the  differential  number 
flux 


Wkd4 

dt 


=  r  kwk 


(ii) 


*.  /  \  _ Pa  0.5  /EminV7"  /«-j\ 

*o(£) — ~~z  :  ^avmin  (  )  (13) 

( pa  0.75).  One-dimensional  (magnetically  field-aligned) 
numerical  modeling  of  electron  acceleration  by  strong 
Langmuir  turbulence  yields  $a(e||)  ~  ejj"1,  consistent  with 
the  ID  scaling  law  for  Langmuir  collapse  [Galeev  et  al., 
1983;  Wang  et  al.,  1997],  A  flat  distribution  of  accelerated 
electrons  is  consistent  with  the  observations  of  Carlson  et 
al.  [1982], 

[40]  The  minimum  energy  emjn  and  the  density  of  par¬ 
ticles  in  the  high-energy  tail  na  are  determined  by  the 
boundary  condition 


^a(emin)  —  4>o(£min) 


(14) 


Figure  8.  Differential  number  fluxes  of  photoelectrons 
(solid  lines)  and  accelerated  electrons  (dash-dotted  lines)  in 
cm_2s_1ster_1eV_1  for  ns  =  1  (line  1)  and  10  (line  2)  cm-3. 


and  the  wave  energy  flux  transferred  by  collapsing  cavitons 
and  absorbed  by  accelerated  particles 


( "ll  \El\  \ 

p  \mi  24ir»07'e  J 


-I*. 


(15) 


Here  Ekm  is  the  wave  amplitude  in  a  cavity  of  the  scale 
length  lmi„  ~  k~'  ~  vmin/wp  when  collapse  is  arrested  and 
$0  is  the  differential  number  flux  of  ionospheric  electrons. 

[41]  For  a  Maxwellian  electron  distribution,  one  has  e^n 
~  20fe.  When  photoelectrons  are  present,  4>0(e)  — »  ^(e)  at 
e  »  Te,  and  emin  is  defined  by 

4h~30 {nJ'lW,)11*  eV  (16) 

valid  for  the  energy  density  of  Langmuir  turbulence  W/  » 
noTpnJm,  [Mishin  and  Telegin,  1986].  Given  ns  =  10  cm-3, 
ne  =  6  103  cm-3,  and  W )  ~  1 0  AneTe,  from  (16)  one  gets 
Emm  —  15  eV.  The  number  of  accelerated  electrons  exceeds 
the  photoelectron  background  in  the  energy  range  e  >  Emm 
by  a  factor  of  ~  10.  For  energies  e  <  Emin,  the 
distribution  is  unaffected.  Figure  8  shows  differential 
number  fluxes  of  photoelectrons  (7)  and  accelerated 
electrons  (13)  for  ns  -  1  and  10  cm-3,  ne  ~  6  105  cm-3, 
and  W,  ~  10_4«ere. 

4.  Discussion 

[42]  Figure  9  shows  the  components  of  the  equilibrium 
volume  emission  rates  ti>f(l)  for  the  period  near  0430  UT. 
To  calculate  the  contributions  of  photoelectrons  (7)  and 
accelerated  electrons  (13),  ns  =  10  cm-3  and  Wt  ~  10~4neT2 
were  chosen,  yielding  Emin  —  15  eV  for  /0  =  7.8  MHz.  As 
expected  from  (4),  the  dissociative  recombination  contribu¬ 
tion  to  the  red  line  is  reduced  much  more  than  that  for  the 
green  line.  Apparently,  heated  electrons  dominate  the  red- 
line  emission,  while  accelerated  electrons  are  the  major 
contributor  to  the  green  line.  Note  a  factor  of  about  2 
difference  between  calculations  with  Fw{7’1)  and  Fth(T{) 
consistent  with  Mishin  et  al.  [2000]  (it  is  ~3.5  for  Te  =  T2). 

[43]  In  the  heated  volume  the  densities  of  excited  oxygen 
atoms  and  of  02  ions  grow.  During  the  onset  period  [O('D)] 
and  [02]  remain  well  below  the  new  equilibrium  values  that 
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Figure  9.  Equilibrium  volume  emission  rate  components 
of  the  red  (a)  and  green  (b)  emissions:  Dissociative 
recombination  with  [O2]o  from  the  background  Te  =  7o 
(thin  solid  line)  and  heated  Te  =  T2  (bold  solid  line)  regions; 
photoelectrons  (thin  dashed  line)  and  accelerated  electrons 
(bold  dashed  line);  0(lS)  cascading  from  the  background 
(circles)  and  heated  (dots)  regions;  thermal  electrons  with 
Te  ~  Tq  (crosses),  7)  (filled  triangles),  and  T2  (open 
triangles);  Maxwellian  electrons  with  Te  =  7)  (asterisks). 


are  of  the  order  of  Q,JLr  and  [02]/a(T2),  respectively. 
Corresponding  time-variations  of  the  contributing  compo¬ 
nents  to  the  total  volume  emission  rate  during  the  onset  can 
be  evaluated  assuming  that  Ajy  <  C,  <  min[l/(«ea2),l/Lr]  — 
15  s  and  Lr  ~  const(f)  at  200-300  km.  Here  Ci  =  t  —  t0  and 
the  heating  time  t0  ~  1  s.  As  a  result,  one  gets 

6[02+  ]  =  [O2+]/[O2+]0-l  a  ».[«(Jo)  ~  c.(7i)]  •  C  (17) 

This  yields 

■HgCO  -  +  tj“  ne[a(T0)  -  «(7'2)]  •  C 

(18) 

n-(0  -  (KqLr  +  0.5r)fR[I,.  +  nea(72)]6[C»^])C 

[44]  Using  ( 1 8)  and  Figure  9,  the  variation  of  the  red-  and 
green-line  emissions  during  the  onset  period  are  calculated 


from  (2).  The  results  are  shown  in  Figure  10  together  with 
the  observed  intensities  (section  2).  There  is  good  apparent 
agreement  with  the  HAARP  observations. 

[45]  We  emphasize  that  the  onset  values  of  the  HF- 
induced  airglow  amount  to  £  1 0%  of  the  saturated  values. 
As  the  reflection  heights  of  7.8  and  6.8  MHz  waves  at 
magnetic  zenith  are  below  the  upper  hybrid  resonance, 
similar  to  Figure  1,  striations  are  not  generated.  Thus  self- 
focusing  due  to  striations  seems  unlikely.  This  indicates  that 
the  ionosphere’s  parameters  in  the  heated  volume  are 
significantly  modified.  Apparently,  a  complete  explanation 
of  artificially  induced  airglow  involves  a  complex  chain  of 
plasma  and  chemical  interactions  within  the  heated  volume 
and  has  yet  to  be  worked  out. 

5.  Conclusion 

[46]  Observations  of  airglow  at  630  nm  (red  line)  and 
557.7  nm  (green  line)  acquired  during  the  February  2002 
campaign  at  the  HAARP  heating  facility  have  been  ana¬ 
lyzed.  During  injections  toward  magnetic  zenith  the  inten¬ 
sities  of  green  and  red  line  emissions  gain  ~5-10  R  within 
~1  s  and  ~20  R  within  ~10  s,  respectively.  We  call  this 
period  the  onset  of  HF-induced  airglow  and  develop  an 
explanatory  model  for  the  growth  of  emissions  from  mag¬ 
netic  zenith.  The  model  accounts  for  the  roles  played  by 
ambient  photoelectrons  and  dissociative  recombination  of 
O*  ions  and  shows  that  heating  and  acceleration  of  back¬ 
ground  electrons,  possibly  by  strong  Langmuir  turbulence, 
dominate  the  airglow  onset. 

[47]  We  also  propose  a  new  scenario  for  the  generation  of 
strong  Langmuir  turbulence  for  injections  outside  the  Spitze 
region,  including  magnetic  zenith,  with  upper  hybrid  waves 
as  the  primary  source  of  Langmuir  waves.  Those  waves 
accumulate  in  the  region  of  small  wavevectors  due  to 
induced  scattering  by  ions  (Langmuir  condensate).  Waves 
in  the  condensate  are  subject  to  the  modulational  instability 
and  collapse  thus  leading  to  strong  turbulence.  This  scenario 


0  5  10  15  20 

time,  s 

Figure  10.  Calculated  intensities  of  the  HF-induced  red 
(solid  line)  and  green  (dashed  line)  emissions  during  the 
onset  period.  The  observed  intensities  are  shown  by  open 
circles  (red)  and  triangles  (green),  respectively. 
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applies  not  only  at  the  reflection  layer  but  also  well  below  it, 
wherever  the  matching  conditions  for  the  parametric  decay 
o  — >  uh  +  Ih  are  met  and  local  values  of  xuh  <  1.  For 
injections  toward  magnetic  zenith  and  heating  frequencies 
>5.4  MHz,  this  occurs  below  the  altitude  of  the  upper 
hybrid  resonance.  The  large  altitude  extent  of  the  strong 
turbulence  region  appears  to  be  the  key  factor  in  exciting 
strong  airglow  at  magnetic  zenith. 
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